Dopant induced ignition of helium nanodroplets in intense few-cycle laser pulses 
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We demonstrate ultrafast resonant energy absorption of rare-gas doped He nanodroplets from in- 
tense few-cycle (~10fs) laser pulses. We find that less than 10 dopant atoms "ignite" the droplet 
to generate a non-spherical electronic nanoplasma resulting ultimately in complete ionization and 
disintegration of all atoms, although the pristine He droplet is transparent for the laser intensities 
applied. Our calculations at those intensities reveal that the minimal pulse length required for 
ignition is about 9fs. 



The ionization dynamics of atomic clusters in intense ul- 
trashort laser pulses has been an active area of research in 
recent years at near-infrared (NIR), vacuum-ultraviolet 
(VUV) and soft X-ray wavelengths [IH3]. The energy 
absorption of NIR light by rare-gas clusters is far more 
efficient as compared to atomic jets or planar solid targets 
when irradiated with similar laser pulses. The underly- 
ing mechanism is resonance absorption by an electronic 
nanoplasma which forms through illumination with the 
light [4-6J. The resonance occurs when the background 
ions of total charge Q occupy a large enough sphere 
of radius R such that the nanoplasma eigenfrequency 
ft = y/Q/R 3 matches the frequency uj of the laser light 
[lj. The resonance condition Q = uj is achieved on the 
time scale of atomic motion (sub- or few-picoseconds) 
since the cluster, i.e. the ions, must expand in order for 
the ions to become sufficiently dilute. This has been ver- 
ified in experiments where pump-probe delay |7J or pulse 
durations are varied |8J, as well as in calculations |5J. 

The fact that this resonance absorption at NIR fre- 
quencies relies predominantly on atomic expansion has 
two consequences: Firstly, it requires time scales of ionic 
motion in spite of its electronic nature. Secondly, it does 
not occur in isolated atoms or those inside the bulk in 
condensed phase. One may ask if there is another possi- 
bility for this resonant coupling of energy from the light 
pulse to the electrons which is not only extremely efficient 
but also fast, bypassing atomic expansion. Indeed this 
can be achieved as theoretically proposed and demon- 
strated |9J, if one uses a two component system with 
a few seed atoms with ionization potentials lower than 
the majority of atoms which belong to the other com- 
ponent. If the spherical symmetry of the system is suit- 
ably broken (e.g. Xe atoms in the center of a He droplet 
irradiated by linearly polarized light) a "cigar shaped" 
nanoplasma is formed. It has two eigenfrequencies which 
bracket the eigenfrequency of a spherical plasma. The 



lower one (along the linear laser polarization) is reso- 
nant with 790 nm laser light right away at typical atomic 
cluster densities leading to almost immediate resonance 
absorption without the need for atomic expansion. This 
allows not only to realize the efficient energy transfer very 
quickly, now limited essentially by the rise time of the 
laser pulse, but also opens the resonant light absorption 
to other forms of matter which do not (Coulomb) ex- 
plode. 

Here, we report on experiments which demonstrate 
this purely electronic resonance absorption with rare-gas 
doped helium droplets illuminated by a few-cycle pulse 
of a duration as short as 10 fs at 790 nm. This rules out 
the influence of any kind of atomic motion. Reducing 
the pulse length even further must eventually prevent 
even this electronic resonance absorption. Indeed, as 
we will demonstrate with calculations, there is a clear 
transition from dopant induced resonance absorption to 
a regime dominated by static field ionization [10] if the 
pulse length falls below a critical value which depends of 
course on the intensity of the laser pulse. 

He nanodroplets can be doped with a well-controlled 
number of rare-gas atoms which aggregate at the center 
of the droplet [TTj fT^ ] . These are exceptionally suitable to 
demonstrate the role of seed atoms in resonance absorp- 
tion, since the laser intensity can be chosen such that the 
light couples almost exclusively to the dopants whereas 
the pristine He droplet is transparent. Interestingly this 
applies to NIR [9j[T3] as well as to X-ray [10] frequencies. 
Despite its transparency, the droplet becomes highly ac- 
tive once the core of seed atoms is ionized. 

In the experiment a beam of He nanodroplets is pro- 
duced by expanding pressurized 4 He gas (70-90 bar) 
through a nozzle (5/im in diameter) maintained at a 
temperature of 15-25 K. By varying the nozzle temper- 
ature in this range the mean number of He atoms per 
droplet is adjusted in the range 10 3 -10 5 . Owing to their 



low internal temperature ~ 370 mK, these droplets are 
in a superfluid state [11] [12]. They are doped by passing 
the skimmed beam of pure droplets through a cylindri- 
cal doping cell that is 3 cm long with collinear apertures 
(0 = 3 mm). The mean number of dopants per nano- 
droplet K is regulated through controlled leaking of the 
desired rare-gas into this cell using a dosing valve (leak 
rate < 10 -10 mbar 1/s). The pressure in the doping cell 
is monitored by a directly attached vacuum gauge. 

Recently, we have simulated the pick-up process of al- 
kali atoms that form high-spin clusters using a Monte- 
Carlo model [14]. Here, we use the same model adapted 
to rare-gas dopants to validate the estimate of doping 
levels and to estimate the loss of He due to evapora- 
tion caused by the doping process: The number K of 
dopant atoms is ascertained from the cell pressure and 
the droplet size according to the modified Poissonian 
pick-up statistics [12J using the semi-empirical formula 
of Kuma et al. [15] which takes evaporation of He atoms 
into account. 

Intense few-cycle laser pulses (~10fs) at a central 
wavelength of 790 nm with peak intensities in the range 
10 14 -10 15 W/cm 2 are generated from a Ti-Sapphire 
based mode-locked laser system. Photo-ions are detected 
by a time-of-flight (TOF) spectrometer in the Wiley- 
McLaren geometry. TOF mass spectra of He ++ ions 
for different mean doping numbers K of Xe in a droplet 
containing 1.5 xlO 4 He atoms are shown in the inset in 
Fig. [11(b). The He ++ peak, which is characteristic for 
ionizing doped He droplets and is absent in mass spectra 
of atomic He gas, is split as a result of the finite angular 
acceptance of the spectrometer into two components due 
to ions with high momenta directed towards and away 
from the detector, respectively. The yield of He ++ and 
He + ions is extracted from the TOF mass spectra by 
integrating over the respective mass peaks. 

Fig(l](a) presents the doping dependence of the yields 
of He~^ and He + at a peak intensity of 7 x 10 14 W/cm 2 
of the laser pulse. A gradual increase of doping num- 
ber K = 1 up to 10 leads to a dramatic step-like in- 
crease in the yields of He + and He ++ ions. We refer 
to this nanoplasma formation as dopant-induced ignition 
(DII). The saturation of ion yields, and equivalently the 
build-up of charge in the nanoplasma, occurs for a criti- 
cal doping number K cr of just ~7(±1.4). The ionization 
of the He droplet is not enhanced any further by adding 
more dopants. Based on classical molecular dynamics 
(MD) [5j [6] we have tailored a theoretical description to 
the present problem of doped He droplets (for details of 
the approach see [9j H3]). The result from our calcu- 
lation shows (Fig(T](b)) qualitatively the same behavior, 
although a direct comparison of the charge to yields is dif- 
ficult since the experiment is an average over the droplet 
size and doping distributions as well as over intensities in 
the focal volume. Moreover, in the theoretical result elec- 
tron recombination is not taken into account. Yet, the 
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FIG. 1. (a) Yields of He ++ and He + ions as a function of 
the mean number of Xe dopants in a nanodroplet containing 
1.5x 10 4 He atoms at a peak laser intensity of 7x 10 14 W/cm 2 . 
(b) Numerical calculation of charge per He atom as a function 
of number of Xe atoms in a droplet containing 4000 He atoms. 
All lines are to guide the eye. Inset: Time-of-flight mass 
spectra of He ++ and He + ions for different numbers of dopant 
Xe atoms (K) as shown in the legend. 
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FIG. 2. Xe doping dependence of He ++ yields at a peak 
laser intensity of 1.5 x 10 14 W/cm 2 for droplets with a 
mean number of He atoms per droplet Nn e =4500 (red) and 
Nne= 15000 (blue). Horizontal bars present the error in esti- 
mating the number of doped Xe atoms. The dash-dot curves 
present the fraction of droplets lost due to complete evapora- 
tion, which is significant only for number of Xe dopant atoms 
K > 20 (grey shading). 



curves agree well qualitatively and even quantitatively 
regarding the critical dopant number at which the sat- 
uration of charging occurs. This agreement points to a 
robust underlying mechanism. 

The critical doping number is remarkably independent 
of the mean droplet size as shown in Fig. [2] having the 
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FIG. 3. The He ++ yields for a droplet containing 15000 He 
atoms as function of the number of Xe dopant atoms (a) for 
different laser intensities I and (b) for different dopant species 
at 7 = 7 x 10 14 W/cm 2 . 



value of 10(±2) for mean droplet sizes of 4500 and 15000 
atoms at a peak laser intensity of 1.5 x 10 14 W/cm 2 . The 
error in estimating the number of dopants using Kuma's 
approach [15 J is also shown. The decrease of doubly 
charged He for increasing K in Figs.[IfJ3] is related to 
the doping process: When dopant atoms are picked up 
sequentially and aggregate inside the droplet, collisional 
as well as binding energy is released into the droplet. 
This leads to an increasing fraction of He droplets that 
completely evaporate at high doping pressure. The cor- 
responding estimate from the Monte- Carlo simulation for 
Xe doping is indicated in Fig. [2] as dashed lines. Clearly, 
the onset of significant nanodroplet destruction in this 
case of Xe doping occurs well beyond the K cr (~ 10) ob- 
served in the experiment. Thus, the process of dopant- 
induced ignition is not influenced by evaporation effects. 

The insensitivity of DII with respect to the droplet 
sizes results simply from the fact that DII sets in with 
the formation of a cigar shaped plasma which is much 
smaller than the droplet itself |9J. In contrast, DII should 
depend on the laser intensity and similarly on the ioniza- 
tion potential of the seed atoms, as well as on the pulse 
length. The former dependences are illustrated with ex- 
perimental results in Fig. [3] As expected, K cr is smaller 
for higher intensities and/or lower ionization potentials. 

We determined numerically the average He charge in a 
droplet Y(K) as a function of the number of Xe dopants 
K. One can see in Fig. [4] that at the intensity of 7 x 
10 14 W/cm 2 a qualitative change of the average charge 
per He atom Y(K) occurs between 7 and 9 fs pulse length: 
For longer pulses DII (solid lines) is operative with a fast 
rise followed by saturation at full ionization of the He 
atoms (Y = 2). For shorter pulse lengths the average 
charge shows a gradual rise Y oc K 3 ! 2 (dashed lines) 
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FIG. 4. Average charge per He atom Y(K) as a function of 
the number of doped Xe atoms for a droplet with 4000 He 
atoms and intensity as in Fig. [I] from MD calculations. The 
different curves distinguish from bottom to top laser pulse 
lengths (half width half maximum) from 3 to 15fs in steps 
of 2fs. Additionally, Y(K) for 10 fs (crossed boxes) is shown, 
see also Fig. [I] The lines are fits of the data points according 
to Eq. |2l (solid) and with the power law Y oc K 3 ^ 2 corre- 
sponding to ionization by a static field for T < 9 fs (dashed) . 
For details, see text. 



which is characteristic for ionization due to the static 
electric field of the laser ionized Xe ions in the center of 
the droplet. 

To derive this dependence we assume that all Xe ions 
constitute a point charge at the center of the droplet and 
that the density of He atoms is uniform within the spher- 
ical droplet. Then the problem is of purely radial nature 
and a sphere exists which contains only doubly ionized 
helium atoms. Its radius is given by i? ++ = y/ qx. e K / E\> S 2 
where qx e is the average charge of each of the K dopant 
atoms and E^^ is the electric field necessary to ion- 
ize the second electron in He through barrier suppres- 
sion (Bethe rule) [TO} ITS] . The number of He atoms 
inside this sphere is given by N ++ = 47rpHe^++/3 = 
47rpHe(#Xe^/^bs2) 3//2 /3, where pn e is the density of He 
in the droplet. We find the number of singly ionized he- 
lium atoms iV_|_ in a similar way and obtain the average 
ion charge per He atom using Y = (iV+ + 2iV++)/iVHej 
where A^e is the total number of He atoms in the droplet. 
This yields 



Y(K) = |£ K/ 2 + iC/ 2 ) &'\ (1) 

where £?bsi , ^bs2 are the first and the second barrier sup- 
pression fields for He. Assuming that qx e does not de- 
pend strongly on the number of Xe atoms K (this is ful- 
filled very well for K between 2 and 15 as seen from the 
MD simulations), we conclude that Y oc K 3 / 2 , i.e. the 
static field ionization due to positively charged dopants 
follows a characteristic power law. 
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FIG. 5. The number of dopants for which the average charge 
per He atom Y in Fig. El has a turning point, Y"(K CV ) — 0. 
The shaded area distinguishes the regime of dopant induced 
ignition (DII) from to the field ionization regime (not shaded), 
at a peak intensity of 7 x 10 14 W/cm 2 . 



demonstrates DII for a pulse length of 10 fs. 

In summary, we have experimentally demonstrated 
that it is possible to transfer energy resonantly from a 
790 nm intense few-cycle pulse to bound electrons with- 
out the need of atomic expansion of the target. This 
brings down the time scale of this extremely efficient but 
relatively slow process from the sub-picosecond regime 
down to a few femtoseconds and, at the same time, allows 
it to be applied to any form of matter that can be suitably 
doped with seed atoms. Finally, the dopant induced igni- 
tion demonstrated here, may also explain the surprising 
enhancement of light absorption in water doped argon 
clusters [T7] , 

Stimulating discussions with J. Tiggesbaumker and 
Th. Fennel, and the valuable assistance of O. Biiner- 
mann in the Monte- Carlo calculations are gratefully ac- 
knowledged. This work is supported by the Deutsche 
For schungsgemeinschaft . 



Although it is clear from Fig.[4]that for increasing pulse 
lengths the average charge per He atom Y{K) is no longer 
well described by the field ionization power law (dashed 
lines), a well defined transition from field ionization to 
DII for increasing pulse length is difficult to extract. To 
determine the minimal pulse length T for which DII dom- 
inates, we introduce a function for the shape of Y(K), 
which can describe both limits, 



Y{K) = Y a 



{K/af 2 - 1 
(K/a) 3 / 2 + ■ 



(2) 



For /3 ^> K/a, Eq. [2] describes the field ionization behav- 
ior, Y(K) oc if 3 / 2 , while for j3 <C K/a, the typical DII 
shape emerges with a sharp onset followed by saturation 
at Yqq. The latter is characterized by a negative second 
derivative Y"(K) while for the field ionization power law 
Y"{K) > holds for all K. Hence, Y"(K CT ) = can 
be interpreted as the conditional equation for the critical 
number of dopant atoms where field ionization dominated 
absorption goes over into DII dominated absorption. Us- 
ing Eq.[2j we get K cr = a(/3/5) 2 / 3 . The parameters (a, (3) 
are obtained for each pulse length T by fitting Eq. [2] to 
the numerically generated data in Fig. [4j K cr exhibits 
an exponential dependence on the pulse length T, as can 
be seen from Fig. [5] The shaded area indicates DII with 
Y" < 0, while below the line of K cr field ionization rules. 
Since physically, at least one dopant atom is necessary 
for ionization, DII dominates for T larger than about 
9fs. This is consistent with the qualitative conclusion 
from Fig. [4] as well as with our experimental result which 
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